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SUMMARY 


Transient  Induced  aerodynamic  data  has  been  analyzed  for  a  subsonic 
VSTOL  aircraft  model  operating  In  a  moving  deck  environment.  Only  the 
single  axis  deck  motions  of  heave,  pitch,  and  roll  have  been  considered. 

The  equations  governing  this  analysis  are  presented  along  with  the 
methods  used  to  reduce  and  analyze  the  data.  Dynamic  effects  are  discussed 
for  each  deck  motion  and  time  history  comparisons  shown  using  developed 
prediction  methodology.  Conclusions  are  presented  for  this  analysis 
along  with  recommendations  for  further  Investigations  into  this  area. 

Appendix  A  presents  a  detailed  discussion  of  the  governing  equations 
along  with  a  description  of  the  computer  programs  utilized.  The  transient 
induced  aerodynamic  data  which  were  analyzed  for  heave,  pitch,  and  roll 
motions  are  presented  in  Appendices  B,  C,  and  D  respectively. 
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INTRODUCTION 


Operating  a  V/STOL  aircraft  in  ground  effect  can  greatly  alter  the 
aerodynamic  forces  and  moments  acting  on  the  aircraft.  By  introducing  a 
moving  ground  plane,  such  as  the  deck  of  a  DD963  class  destroyer,  the 
complexity  of  the  resulting  flow  field  is  significantly  increased.  This 
added  complexity  makes  it  even  more  difficult  to  perform  successful 
shipboard  takeoff  and  recovery  operations.  In  order  to  fully  utilize 
the  potential  shipboard  operational  capabilities  of  V/STOL  aircraft  a 
better  understanding  of  these  phenomena  is  required. 

To  accomplish  this,  the  effects  of  basic  single  axis  deck  motions 
must  be  examined.  Analysis  of  these  independent  deck  motions  (heave, 
pitch,  roll)  enable  the  development  of  prediction  technqiues  applicable 
to  a  preliminary  design  environment.  These  methods  provide  the  designer 
a  tool  for  estimating  the  transient  aduced  aerodynamics  acting  on  a 
V/STOL  aircraft  in  a  moving  deck  environment. 

Data  for  the  three-fan  subsonic  V/STOL  configuration  of—ref erenee  (a) 
was  selected  for  study. three  view  drawing  of  the  full  scale  configu¬ 
ration  is  shown  in  Figure""!*^  As  evident  in  this  figure,  the  aircraft 
has  a  low  wing  and  high  tail  with  two  lift /cruise  fans  mounted  above  the 
wings  and  one  lift  fan  in  the  forward  fuselage.  The  method  of  analyzing 
this  data  in  terms  of  single  axis  deck  motions  will  be  discussed  in 
detail  in  the  following  sections. 


APPROACH 


This  analysis  was  limited  to  sinusoidal  deck  motions  of  heave, 
pitch,  and  roll  only.  Combined  motions  such  as  heave  and  roll  were  not 
considered  due  to  the  complexity  that  they  introduce.  These  combined 
motions,  however,  should  be  considered  in  future  analyses  of  transient 
induced  aerodynamics  acting  on  a  VSTOL  aircraft. 


The  dynamic  effects  caused  by  a  moving  deck  were  investigated 
through  examination  of  the  induced  forces  or  moments  in  terms  of  frequency, 
amplitude,  and  height-above-deck.  For  each  deck  motion  the  induced 
force  or  moment  was  assumed  to  be  a  function  of  deck  position  and  rate. 

In  the  cases  of  pitch  and  roll  the  effect  of  height- above-deck  was 
also  included  as  shown  below: 
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The  equations  defining  the  deck  positions  and  attitudes  at  any  time,  which  were 
differentiated  to  obtain  the  desired  rates,  are  presented  below: 


h  =  hQ 

+  Ah  sin(oit) 

2 

(4) 

9  -  0O 

+  A0  sin  (wt) 

2 

(5) 

<t>  =  $o 

+  A£_  sin  (ait) 

(6) 

2 


Figure  2  defines  the  various  parameters  of  these  equations.  The  deck 
position  and  rate  are  used  in  a  correlative  analysis  technique  to  determine  the 
transient  induced  effects  as  illustrated  in  Figure  3  ((a)  through  (e)).  This 
example  is  limited  to  the  transient  induced  lift  resulting  from  variations  in 
deck  heave  only;  this  being  typical  of  the  method  applied  to  the  cases  for  pitch 
and  roll. 

Through  selection  of  a  deck  rate  (h/de)  the  corresponding  position  (h/de) 
and  time  (t^,  t2,  •••)  were  obtained.  Figure  3(a).  Knowing  the  measured  induced 
lift  as  a  function  of  time,  Figure  3(b),  and  the  corresponding  deck  position, 
the  induced  lift  was  cross-plotted  with  position.  Since  the  deck  motion  was 
sinusoidal  these  data  occurred  in  pairs  about  the  neutral  point  (ho/de)  for  the 
specified  deck  rate.  Figure  3(c).  By  varying  deck  rate,  an  envelope  of  induced 
lift  versus  position  was  developed  corresponding  to  a  given  set  of  test  para¬ 
meters,  i.e.  ,  deck  neutral  point,  ho/de,  frequency,  ui,  and  amplitude,  Ah/de. 

From  this  map  a  quasi-static  (h/de  *  0)  induced  lift  was  obtained  in  order  to 
predict  the  actual  transient  induced  effects.  Figure  3(d).  Comparisons  of  these 
predicted  transient  effects  with  the  actual  test  data.  Figure  3(e),  are  pre¬ 
sented  in  subsequent  figures  for  selected  test  conditions. 

RESULTS 

Each  of  the  three  deck  motions  will  be  discussed  in  terms  of  the  accompany¬ 
ing  static  and  dynamic  aerodynamic  characteristics.  Comparisons  of  actual  time 
histories  of  dynamic  data  with  predicted  time  histories  will  also  be  presented. 

HEAVE 


For  the  heaving  deck  motion,  induced  lift  is  the  parameter  which  is 
examined.  Figures  4,  5,  and  6  show  induced  lift  plotted  versus  height  above* 
deck  for  frequencies  of  1,  2,  and  3  Hz.  As  noted  on  these  figures,  the  open 
symbols  correspond  to  positive  rates  (model  heaving  away  from  deck)  while  the 
closed  symbols  represent  negative  rates  (model  heaving  toward  the  deck).  For 
clarity,  only  the  data  run  at  the  largest  amplitude  is  plotted  for  each  fre¬ 
quency  although  all  amplitudes  tested  were  considered  in  the  analysis.  Also, 
the  neutral  points  for  these  figures  correspond  to  the  data  run  presented  and 
not  necessarily  to  the  other  amplitudes. 


DECK  HEAVE  OSCILLATION _  |  MEASURED  INDUCED 


Typical  Correlaiwe  Analysis  Procedure 
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In  Figure  4,  ou  =  1Hz,  the  dynamic  data  differs  from  the  static  data 
for  both  positive  and  negative  rates.  This  is  also  the  case  for  w  *  2Hz 
and  u>  =  3Hz  as  shown  in  Figure  5  and  6  respectively.  In  each  case  the 
dynamic  data  shows  a  hysteresis  type  of  effect  depending  on  the  rate 
being  positive  or  negative.  For  negative  rates,  the  dynamic  results 
show  an  increase  in  the  static  data  over  all  values  of  h/de  for  u  ■  1 
and  to  =  3Hz.  As  shown  for  w  =  2Hz  for  negative  rates  the  dynamic  results 
indicate  a  decrease  in  the  static  AL/T  for  h/de  less  than  1.1  and  an 
increase  for  h/de  greater  than  1.1.  In  each  case  for  positive  rates  the 
results  indicate  a  decrease  in  the  static  AL/T  up  to  a  value  of  h/de 
between  1.25  and  1.50.  After  this,  an  increase  is  seen  up  to  an  h/de  of 
3.5  which  was  the  maximum  value  tested.  This  increase  gets  smaller  as 
the  frequency  increases  from  oi  =  1  to  3  Hz  with  values  close  to  static 
value  at  u  =  3 Hz.  The  hysteresis  effect  described  is  consistent  with 
the  findings  of  reference  (a)  and  is  attributed  to  the  compression 
effect  associated  with  the  fountain.  In  the  case  of  a  positive  rate, 
where  the  deck  is  heaving  away  from  the  model,  a  larger  value  of  AL/T  is 
evident  as  compared  to  a  negative  deck  rate  with  the  same  magnitude. 

This  indicates  that  positive  rates  amplify  the  suckdown  while  negative 
rates  amplify  the  fountain  force. 

Data  from  Figures  4,  5,  and  6  along  with  data  for  the  amplitudes 
which  are  not  plotted  were  used  to  develop  a  curve  of  induced  lift 
versus  height-over-deck  for  a  zero  value  of  deck  rate.  This  curve  is 
shown  as  Figure  7  along  with  the  static  test  data.  As  seen  in  this 
figure,  from  an  h/de  of  0.5  to  1.0  the  curve  indicates  a  decrease  in 
AL/T  as  compared  to  the  static  data.  Above  h/de  =  1.0,  the  curve  shows 
an  increase  in  AL/T  for  all  values  up  to  the  maximum  of  h/de  =3.5. 

This  curve  of  zero  deck  rate  was  used  to  estimate  the  induced  lift 
as  a  function  of  h/de.  By  comparing  these  estimates  with  actual  time 
history  data  the  dynamic  effects  could  be  examined.  These  time  histories 
are  shown  in  Figures  8a  through  8h  with  the  neutral  point,  amplitude, 
and  frequency  Indicated.  Also  shown  are  the  estimated  values  from 
Figure  7  as  indicated  by  the  dotted  lines.  Figures  8a  through  8d 
(co  =  2Hz)  are  presented  first  since  a  more  complete  data  base  was  avail¬ 
able  for  this  frequency.  Very  good  agreement  is  seen  for  an  h0/de  =0.8 
over  the  time  segment  shown  with  peak  values  matching  very  well  and  no 
apparent  time  shift. 

For  hQ/de  =  2.0,  Figures  8b  through  8d,  the  predicted  values  are 
within  1.2%  of  the  actual  data.  As  the  amplitude.  Ah/de,  is  increased 
the  response  shows  a  larger  suckdown  effect.  For  Ah/de  *  1.0  the  smallest 
value  of  AL/T  is  1.0%  while  at  Ah/de  »  3.0  a  value  of  -6.0 %  is  obtained. 
Also,  the  response  becomes  smoother  and  more  sinusoidal  as  Ah/de  is 
increased. 

Figures  8e  and  8f  show  time  histories  for  w  =  1Hz  and  3Hz  at  an 
ho/de  =0.8  and  Ah/de  =0.6  respectively.  In  Figure  8e  a  good  correlation 
is  shown  which  is  within  1.0%  of  the  actual  values  for  the  time  segment 
presented.  It  should  be  noted  that  these  time  histories  are  representative 
responses  form  the  original  5  second  data  runs  of  reference  (a).  For 


PREDICTED 


Figure  8b.  Induced  Lift  Time  History  Comparison  (Run  90.5) 


Figure  8d.  Induced  Lift  Time  History  Comparison  (Run  90.2) 


PREDICTED 


Figure  8f.  Induced  Lift  Time  History  Comparison  (Run  89.3) 


NADC-79260-60 


to  =  3Hz,  figure  8f,  good  agreement  is  seen  with  the  actual  data  where 
peak  values  are  within  1.0%.  For  increasing  values  of  AL/T  and  time  the 
predicted  values  are  within  1.0%  while  a  maximum  difference  of  2.0%  is 
seen  for  decreasing  values  of  AL/T  and  increasing  time. 

As  indicated  by  figures  8e  and  8f  for  to  =  1  and  3Hz  respectively 
and  figure  8a  for  to  *  2Hz,  frequency  does  not  appear  to  effect  induced 
lift  at  a  neutral  point  of  ho/de  *  0.8.  In  each  case  the  maximum  positive 
value  of  AL/T  is  approximately  1.4%  while  maximum  negative  values  are 
about  -6.0%  AL/T. 

A  slight  phase  shift  is  evident  for  frequencies  of  to  ■  1  and  3Hz  although 
peak  values  show  fair  correlation  as  indicated  in  figures  8g  and  8h. 

These  data  correspond  to  a  neutral  point  of  ho/de  =2.0  and  an  amplitude 
of  Ah/de  =  3.0  as  does  the  data  of  figure  8d  for  a  frequency  of  to  =  2Hz. 
However,  no  phase  shift  is  evident  for  this  frequency  (to  =  2Hz)  or  any 
of  the  other  data  runs  except  the  two  previously  mentioned  cases  (figures  8g 
and  8h).  This  indicates  that  the  phaseshift  probably  resulted  from  the 
inability  to  precisely  define  the  initial  time  used  for  plotting  these 
data  runs. 

Several  observations  can  be  summarized  with  regard  to  the  heaving 
motion  results  described  above.  First,  a  hysteresis  effect  is  evident 
dependent  upon  the  direction  of  the  motion.  For  a  positive  deck  rate,  a 
larger  value  of  AL/T  is  evident  as  compared  to  a  negative  rate  with  the 
same  magnitude.  Second,  developing  a  curve  of  induced  lift  versus 
height-above-deck  for  a  zero  deck  rate  to  estimate  the  dynamic 
response  correlates  reasonably  well.  Third,  increasing  the  amplitude 
from  1.0  to  3.0  showed  an  increase  in  suckdown  of  nearly  7.0%  indicating 
that  amplitude  can  have  a  large  effect.  Fourth,  frequency  showed  only  a 
small  increase  in  suckdown  for  a  given  amplitude  and  neutral  point. 

Fifth,  a  slight  phase  shift  was  evident  for  frequencies  of  a  *  1  and  3Hz 
at  a  neutral  point  of  ho/de  =  2.0  and  is  attributed  to  the  difficulty  in 
precisely  identifying  an  initial  time  for  plotting  these  data  runs. 

PITCH 

For  the  pitching  deck  motion,  the  induced  pitching  moment  is  the 
parameter  examined.  Figures  9  through  12  show  induced  pitching  moment 
(AM/Tc)  plotted  versus  pitch  angle  (0)  for  id  =2Hz  at  various  heights 
and  amplitudes.  Positive  and  negative  rates  (0)  are  shown  on  these 
figures  by  open  and  closed  symbols  respectively,  along  with  the  static 
test  data  indicated  by  the  solid  line.  The  effect  of  frequency  could 
not  be  determined  from  the  available  data  since  only  u»  =  2Hz  was  tested. 

At  an  h/de  =  0.8,  figure  9  shows  significantly  larger  negative 
values  of  pitching  moment  for  the  dynamic  data  where  the  magnitude  is 
approximately  twice  the  static  value.  For  A0  =  4°  a  small  hysteresis 
effect  is  present  while  A0  =  12®  does  not  show  this  effect.  The  most 
pronounced  hysteresis  effect  is  shown  for  the  case  of  A0  »  20°  where 
differences  between  positive  and  negative  rates  are  as  much  as  2.0%.  It 
is  interesting  to  note  that  for  negative  pitch  angles  *-he  positive  rates 
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Figure  8h.  Induced  Lift  Time  History  Comparison  (Run  90.3) 
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show  a  more  positive  value  while  the  reverse  is  true  for  positive  pitch 
angles.  Increasing  amplitude  does  not  show  a  significant  effect  on  the 
pitching  moment  for  this  h/de  as  the  data  falls  within  a  small  band  for 
the  three  amplitudes  shown. 

Figure  10  shows  the  induced  pitching  moment  at  h/de  =  2.0  for  the 
three  amplitudes  of  4®,  12°,  and  20®.  As  expected,  the  hysteresis 
effect  is  very  pronounced  since  the  fountain  is  the  strongest  at  this 
height.  At  A©  *  4®  positive  rates  show  more  positive  values  over  all 
pitch  angles  than  do  the  negative  rates.  The  reverse  is  true  for  amplitudes 
of  12®  and  20®  with  differences  as  much  as  4.0  AM/Tc  evident.  For  this 
height  the  static  data  shows  a  significant  increase  over  the  values  at 
h/de  =  0.8  where  all  values  shown  are  negative.  The  dynamic  data  also 
shows  an  increase  on  the  order  of  3.0%  as  the  height  is  increased  from 
0.8  to  2.0  as  indicated  in  figures  9  and  10.  It  should  also  be  pointed 
out  that  for  an  h/de  =2.0  the  static  values  become  positive  over  almost 
the  entire  range  of  pitch  angles  shown. 

As  indicated  in  figure  11,  increasing  h/de  to  5.0  produces  significant 
changes  in  pitching  moment.  The  static  data  is  positive  for  both  positive 
and  negative  pitch  angles  while  the  dynamic  data  shows  an  unexpected 
increase  at  A©  =  20®.  This  is  also  the  case  for  an  h/de  =  8.0  as  indicated 
in  figure  12.  At  heights  of  5.0  and  8.0  the  induced  pitching  moment 
should  be  very  small  in  magnitude  and  close  to  the  static  values  since 
the  model  is  essentially  out  of  ground  effect.  For  A©  =  4°  and  12° 
these  data  show  an  almost  constant  value  of  -0. 025  for  induced  pitching 
moment.  This  fact,  along  with  the  data  for  A©  =  20°  being  so  close  to 
the  static  value  at  both  h/de  =5.0  and  8.0  indicates  a  zero  shift  in 
the  data.  This  shift  may  also  apply  to  figures  9  and  10  where  large 
differences  are  seen  between  the  static  and  dynamic  responses.  In  order 
to  develop  a  curve  for  a  zero  value  of  pitch  rate  to  estimate  the  dynamic 
response  it  is  necessary  to  account  for  this  zero  shift.  To  accomplish 
this  the  response  for  A©  =  20®  at  both  h/de  =5.0  and  h/de  =  8.0  was 
assumed  to  fall  at  a  constant  level  of  -0.025  to  be  consistent  with  the 
data  for  A©  =  4°  and  A©  =  12®.  Also,  the  effect  of  height  must  be 
included  since  significant  differences  in  pitching  moment  are  seen  at 
various  values  of  h/de.  Figure  13  presents  the  curves  developed  for  the 
four  values  of  h/de  which  were  tested.  The  height  effect  was  accounted 
for  in  this  way  to  show  that  at  lower  values  of  h/de  this  effect  is  much 
stronger  than  at  the  higher  values. 

Figure  13  was  used  to  predict  the  induced  pitching  moment  as  a 
function  of  pitch  angle  and  height.  These  predictions  were  compared 
with  actual  time  history  data  and  are  shown  in  figures  14a  through  141 
with  predicted  values  indicated  by  the  dotted  line. 

Figures  14a  through  14c  present  induced  pitching  moment  versus  time 
for  h/de  =  0.8  and  amplitudes  of  4,  12,  and  20®  respectively.  Good 
correlation  is  seen  between  the  actual  and  predicted  values  for  each 
amplitude  with  differences  of  less  than  1.0%  AM/Tc".  As  the  amplitude  is 
increased  the  magnitude  of  the  dynamic  response  shows  larger  fluctuations 
of  induced  pitching  moment. 
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Static  and  Dynamic  Induced  Pitching  Moment,  h/de 


Figure  11.  Static  and  Dynamic  Induced  Pitching  Moment,  h/de 
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Figure  13.  Induced  Pitching  Moment  vs.  Pitch  Angle  for  Zero  Deck  Rat 
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Induced  Pitching  Moment  Time  History  Comparison  (Run  82.2 
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For  h/de  *  2.0,  figures  I4d  through  14f  show  induced  pitching 
moment  versus  time  for  the  three  amplitudes  tested.  A  fair  correlation 
is  shown  between  the  predicted  and  actual  responses  at  each  amplitude. 

At  A9  *  4°  differences  of  up  to  1.5%  are  evident  while  at  A9  *  20°  an 
error  of  2.5%  can  be  seen.  These  figures  also  show  that  increasing 
amplitude  has  the  same  effect  on  induced  pitching  moment  at  this  h/de  as 
it  did  at  h/de  -0.8,  where  larger  fluctuations  in  the  magnitude  of  the 
dynamic  response  are  seen  for  Increasing  amplitudes. 

Induced  pitching  moment  versus  time  for  an  h/de  -  5.0  is  shown  in 
figures  14g  through  14i.  The  dynamic  response  at  this  height  should  be 
nearly  zero  since  the  model  is  essentially  out  of  ground  effect.  Figures  14g 
and  14h  show  good  agreement  between  predicted  and  actual  values  although 
the  magnitudes  are  on  the  order  of  -0.025  AM/T5-.  These  values  are  not 
zero  due  to  the  zero  shift  which  has  been  previously  discussed.  Another 
result  of  this  shift  is  seen  in  figure  14i  where  a  poor  correlation 
exists  between  the  predicted  and  actual  values  even  though  the  actual 
response  is  very  small  as  expected. 

In  figures  14j  through  141,  for  an  h/de  =  8.0,  characteristics 
similar  to  those  just  discussed  for  an  h/de  =5.0  are  present.  A  good 
correlation  is  seen  for  amplitudes  of  4°  and  12°  although  the  magnitudes 
are  again  on  the  order  of  -0.025  AM/Tc\  For  A9  =  20°  (figure  141)  a 
poor  correlation  is  seen  even  though  the  magnitude  of  the  actual  response 
is  nearly  zero  as  anticipated. 

Several  observations  can  now  be  summarized  with  regard  to  the 
pitching  motion  results  described  above.  First,  reasonably  good  correlation 
is  seen  between  the  actual  data  and  the  predicted  values  using  a  curve 
of  AM/Tc  versus  pitch  angle  for  a  zero  pitch  rate.  Second,  a  hysteresis 
effect  is  present  in  the  data  which  is  most  pronounced  at  an  h/de  =2.0 
where  the  fountain  is  the  strongest.  Third,  height  shows  a  large  effect 
on  AM/Tcf  at  the  lower  values  of  h/de,  0.8  and  2.0,  while  at  heights 
greater  than  5.0  little  or  no  effect  is  seen.  Fourth,  frequency  effects 
could  not  be  determined  from  these  data  since  w  =  2Hz  was  the  only 
frequency  tested.  Fifth,  increasing  the  amplitude  for  a  given  height 
produces  larger  fluctuations  in  the  magnitude  of  the  actual  response. 

Sixth,  a  zero  shift  is  apparent  in  the  data  for  h/de  =5.0  and  8.0  which 
causes  large  differences  between  actual  and  predicted  values  at  an 
amplitude  of  A9  =  20° . 

ROLL 


For  the  rolling  deck  motion  the  parameter  examined  is  the  induced 
rolling  moment,  Af/Tb.  Figures  15  and  16  show  induced  rolling  moment 
versus  roll  angle  for  the  various  amplitudes,  frequencies,  and  height- 
above-deck  conditions  which  were  tested.  The  open  symbols  indicate  a 
positive  roll  rate  ($),  while  negative  rates  are  shown  as  closed  symbols. 
The  static  data  for  each  figure  is  indicated  by  the  solid  line. 
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Figure  1 4e .  Induced  Pitching  Moment  Time  History  Comparison  (Run  82.3) 
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Figure  141.  Induced  Pitching  Moment  Time  History  Comparison  (Run  82.4) 
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I 4 I .  Induced  Pitching  Moment  Time  History  Comparison  (Run  85.3) 


NADC-79260-60 


<j 

o> 


o 

"3* 


a 

Cz3 
H 
hJ  o 

<  M 

p  a 
H  W 

CJ  od 

<  CL 


•  OC 

«j  oj 

:  c  -a 

4  CO  *<T 


I  II  II 


4-8 


45 


*3- 

ii 

3C 


Induced  Pitching  Moment  Time  History  Comparison  (Run  87.1) 
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Figure  14k.  Induced  Pitching  Moment  Time  History  Comparison  (Run  87.2) 


15.  Static  and  Dynamic  Induced  Rolling  Moment 


Figure  16.  Static  and  Dynamic  Induced  Rolling  Moment,  h/de  =  2.0 
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Figure  15  shows  Che  variation  of  induced  rolling  moment  with  changes 
in  roll  angle  (A4> )  for  w  =  2Hz  at  an  h/de  =  0.8.  The  dynamic  data 
differs  from  the  static  data  for  each  angular  displacement  at  both 
positive  and  negative  rates  as  shown.  For  $  =  0°  the  static  value  of  Ajf/Tb 
is  zero  while  the  dynamic  value  shows  a  shift  of  about  0.15%.  Both 
static  and  dynamic  data  appear  fairly  symmetric  about  $  =  0°  as  expected 
for  this  motion.  Varying  A$  from  4°  through  20°  does  not  show  a  large 
effect  on  rolling  moment  as  indicated.  The  data  for  A4>  *  48,  12°,  and 
20°  fall  within  a  relatively  small  band  similar  in  shape  to  the  static 
curve  but  shifted  up  0.15%  for  both  positive  and  negative  roll  angles. 
Although  there  are  differences  in  rolling  moment  for  positive  and  negative 
rates,  these  differences  are  small  in  magnitude.  A  small  hysteresis 
effect  is  present  in  this  figure,  but  is  not  a  pronounced  effect  at  this 
h/de.  For  h/de  “2.0  and  A<j>  =  30°,  as  shown  in  figure  16,  the  hysteresis 
effect  is  very  pronounced  at  both  frequencies  indicated.  It  should  be 
noted  that  for  w  *  2Hz  the  negative  rates  indicated  by  the  closed  symbols 
show  a  higher  value  of  induced  rolling  moment  than  do  the  positive 
rates.  This  is  true  except  at  the  extreme  roll  angles  where  a  slight 
reversal  is  indicated.  For  ui  =  3Hz  the  opposite  trend  is  indicated 
where  the  positive  rates  show  a  higher  rolling  moment  than  do  the  negative 
rates.  The  difference  in  magnitude  for  positive  and  negative  rates  at 
this  h/de  is  much  larger  than  for  an  h/de  =  0.8.  This  is  due  to  the 
fact  that  for  this  configuration  the  fountain  is  strongest  at  an  h/de  =2.0 
and  therefore  the  largest  variations  should  be  seen  at  this  height.  The 
results  for  pitching  moment  also  showed  the  largest  variation  at  an 
h/de  =  2.0,  although  the  overall  magnitudes  are  on  the  order  of  +0. 35% 
which  is  very  small. 

By  combining  figures  15  and  16  for  rolling  moment  at  the  various 
amplitudes,  frequencies,  and  heights  indicated,  a  curve  of  induced 
rolling  moment  versus  roll  angle  was  obtained  corresponding  to  a  Zero 
value  of  deck  roll  rate.  Figure  17  shows  this  curve  indicated  by  the 
dotted  line  along  with  the  static  test  data  for  both  h/de  =0.8  and  2.0. 

As  indicated  in  this  figure  the  dotted  curve  does  not  pass  through  the 
origin.  For  this  rolling  motion  a  symmetry  about  the  origin  was  expected 
as  shown  for  the  static  data,  but  a  shift  of  about  0.15%  is  indicated. 

This  shift  was  investigated  and  the  probable  cause  is  a  bias  in  setting 
the  neutral  point  statically  and  dynamically.  Assuming  this  shift  is 
due  to  measurement  error,  the  zero  curve  can  be  lowered  0.15%  to  pass 
through  the  origin  and  compared  to  the  static  data  at  h/de  =0.8  and 
2.0.  This  shifted  zero  curve  was  used  to  predict  the  rolling  moment  as 
a  function  of  roll  angle  and  compared  with  actual  time  history  data. 

These  comparisons  are  shown  in  figures  18a  through  18e  with  predicted 
values  indicated  by  the  dotted  line.  Since  the  values  for  rolling 
moment  are  very  small  the  discussion  of  these  figures  will  focus  on 
trends  and  not  on  percent  differences  between  predicted  and  actual 
values.  It  should  be  noted  that  the  effect  of  frequency  on  the  induced 
rolling  moment  at  this  h/de  cannot  be  determined  since  only  one  frequency 
was  tested. 

Figures  18a  through  18c  show  induced  rolling  moment  versus  time  for 
<*>  *  2Hz ,  h/de  =  0.8,  and  at  amplitudes  of  4®,  12°  and  20°.  The  predicted 
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Figure  18a.  Induced  Rolling  Moment  Time  History  Comparison  (Run  92.1) 
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values  correlate  reasonably  well  with  the  actual  response  for  each 
amplitude.  As  the  amplitude  is  increased  from  4°  to  12°  the  positive 
peak  values  show  a  large  increase  while  negative  peak  values  remain 
fairly  constant.  A  further  increase  from  12°  to  20°  only  shows  an 
increase  in  negative  peak  value. 

Figures  18d  and  18e  show  induced  lift  for  h/de  =  2.0,  A<}>  *  30°,  and 
for  two  frequencies  (w  =  2  and  3Hz).  In  both  cases  a  slight  time  shift 
seems  apparent,  however  this  shift  is  attributed  to  the  difficulty  in 
precisely  identifying  an  initial  time  used  to  plot  the  time  history. 

Peak  values  correspond  well  with  the  actual  data  and  frequency  changes 
show  little  or  no  effect  on  rolling  moment.  Although  the  fountain  is 
the  strongest  at  h/de  =  2.0,  no  appreciable  difference  in  magnitude  is 
seen  between  the  response  at  h/de  =  0.8,  A$  =  20°  and  at  h/de  =  2.0, 

A<t  =  30°.  Based  on  this  limited  data  base,  height  does  not  seem  to  have 
a  large  effect  on  induced  rolling  moment. 

Several  observations  can  be  summarized  in  regard  to  the  rolling 
deck  motion  results  described  above.  First,  the  magnitude  of  the  response 
for  any  of  the  conditions  tested  is  very  small  (less  than  0.5%).  Second, 
the  dynamic  response  showed  a  shift  of  approximately  0.15%  above  the 
static  data  and  is  attributed  to  a  bias  in  setting  the  neutral  point. 

Third,  based  on  the  limited  amount  of  data  available,  frequency  and 
height  do  not  appear  to  have  a  large  effect  on  induced  rolling  moment. 
Fourth,  an  increase  in  amplitude  does  not  have  a  large  effect  for  the 
various  amplitudes,  frequencies,  and  heights  which  were  tested.  Fifth, 
a  slight  time  shift  was  observed  in  several  of  the  figures  and  is  attributed 
to  the  difficulty  in  identifying  precisely  an  initial  time  from  the 
data. 


CONCLUSIONS  AND  RECOMMENDATIONS 


Transient  induced  aerodynamic  data  for  a  3-fan  subsonic  VSTOL  model 
operating  in  a  moving  deck  environment  have  been  analyzed  in  terms  of 
single  axis  deck  motions.  These  include  heaving,  pitching,  and  rolling 
motions  only.  Combined  deck  motions,  such  as  heave  and  pitch,  were  not 
considered  for  this  first  correlation  analysis  due  to  the  added  complexity 
that  they  introduce. 

A  method  has  been  developed  to  predict  the  dynamic  response  of 
induced  forces  and  moments  based  upon  a  zero  value  of  deck  rate  for  each 
motion.  In  general,  good  correlation  is  seen  between  the  actual  and 
predicted  values.  The  following  represent  the  more  pertinent  conclusions 
which  have  been  reached  as  a  result  of  this  analysis. 

HEAVE 


•  A  hystersis  effect  is  present  for  this  motion  where  the  induced 
lift  loss,  AL/T,  is  greater  for  the  deck  heaving  away  from  the  model  as 
compared  to  the  deck  heaving  towards  the  model. 
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•  For  a  given  amplitude  and  neutral  point,  frequency  did  not 
have  a  large  effect  on  the  induced  lift. 

•  Increasing  the  amplitude  increased  significantly  the  suckdovm 
effect.  As  much  as  7.0%  suckdown  was  seen  for  an  increase  of  amplitude 
from  1.0  to  3.0  h/de. 

•  Predictions  of  transient  lift  variations  based  on  a  quasi¬ 
static  (zero  deck  rate)  assumption  correlate  within  about  1.5%  of 
actual  test  results. 

PITCH 

•  Height  has  a  significant  effect  on  induced  pitching  moment  at 
lower  values  of  h/de.  At  values  of  h/de  *  5.0  and  8.0,  height  effects 
become  negligible  because  the  model  is  out  of  ground  effect. 

•  A  hysteresis  is  present  and  most  noticeable  at  h/de  -  2.0 
where  the  fountain  is  the  strongest. 

•  Frequency  effects  could  not  be  determined  because  of  the 
limited  data  available. 

•  Increasing  the  amplitude  of  the  sinusoidal  motion  increases 
the  magnitude  of  the  actual  response  by  up  to  1.5%  AM/Tc\ 

•  For  h/de  =5.0  and  8.0  a  zero  shift  is  apparent  in  the  data 
and  is  on  the  order  of  about  -0.025  AM/Tr. 

ROLL 

•  The  magnitude  of  induced  rolling  moment  for  any  of  the  heights, 
amplitudes,  and  frequencies  tested  was  less  than  0.5%. 

•  A  shift  of  approximately  0.15%  was  observed  between  the  dynamic 
response  and  the  static  data  at  both  heights  tested  due  to  a  bias  in 
setting  the  neutral  point  (<£0). 

•  Amplitude  did  not  have  a  large  effect  on  induced  rolling 
moment  for  any  of  the  conditions  tested. 

•  Frequency  did  not  have  a  significant  effect  on  the  induced 
rolling  moment  based  on  the  limited  amount  of  data  available. 

These  conclusions  show  that  a  moving  deck  can  produce  substantial  changes 
in  induced  lift  and  pitching  moment  for  this  configuration.  The  major 
contribution  for  heaving  motion  is  the  amplitude  where  as  much  as  a  7.0% 
difference  in  induced  lift  was  seen  for  two  different  values  of  Ah/de. 

In  the  case  of  pitching  motion,  the  driving  factor  is  height  over  deck 
with  a  3.0%  difference  in  pitching  moment  evident  between  an  h/de  =0.8 
and  2.0.  This  3.0%  difference  can  be  Important  to  a  pilot  in  hover  who 
is  trying  to  maintain  trim. 
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This  analysis  has  shown  that  a  quasi-static  curve,  as  developed 
for  each  deck  motion  described,  can  provide  a  reasonable  estimate  of 
the  dynamic  response.  As  more  data  becomes  available  it  is  recommended 
that  this  approach  be  modified  as  required  to  improve  estimation  accuracy 
at  the  preliminary  design  stage. 

It  is  also  recommended  that  the  data  of  reference  (a)  for  the 
supersonic  conf iguration  be  examined  in  a  similar  manner  to  facilitate  a 
comparison  of  the  subsonic  and  supersonic  configurations.  This  comparison 
would  be  especially  valuable  since  the  supersonic  configuration  is 
dominated  by  suckdown  and  the  subsonic  configuration  is  more  fountain 
dominated. 
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This  appendix  outlines  the  formulation  of  the  governing  equations 
for  sinusoidal  deck  motions.  The  reduction  of  the  reference  (a)  data 
into  a  form  suitable  for  analysis  as  well  as  the  computer  program  developed 
to  analyze  the  data  will  also  be  discussed. 

GOVERNING  EQUATIONS 

In  the  heaving  direction,  induced  lift  was  assumed  to  be  a  function 
of  deck  position  and  rate, 

ft-  -  fx(h)  +  f2(h)  (1) 

Since  the  deck  motion  is  sinusoidal,  the  position  of  the  aircraft 
at  any  time  (t)  can  be  obtained  by  the  neutral  point  plus  a  value  determined 
by  the  sinusoidal  motion.  This  can  be  written  in  equation  fora  as 
follows: 


h  =  h0  +  sin  (ut)  (2) 

The  equivalent  nozzle  jet  exit  diameter,  de,  was  used  to  nondimens ionalize 
equation  (2)  in  order  to  be  consistent  with  the  data  of  reference  (a). 


h_ 

de 


ho 

de 


(ait) 


(2a) 


Using  equation  (2a)  the  deck  rate  was  computed  as  a  function  of  frequency, 
amplitude,  and  neutral  point  which  are  the  parameters  varied  in  the 
reference  (a)  data.  Differentiating  with  respect  to  time  yields. 


h_ 

de 


oi  cos  (nit) 


(3) 


Rewriting  equation  (2a), 


sin  (uit) 


(—  — ' ) 

\  de  -  de  / 


Ah 

de 


(4) 


the  aquation  for  deck  rate  is  obtained  and  can  be  written  in  the  following 
form. 


_h_ 

de 


(5) 
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In  the  pitch  direction,  the  induced  pitching  moment  was  assumed  to 
be  a  function  of  pitch  angle,  pitch  rate,  and  deck  position  as  follows, 

fjf  ”  Si  W)  +  g2(®)  +  83(h)  (6) 

Pitch  angle  and  rate  were  computed  as  follows, 

Afl 

9  ’  90  +  y  sin  (wt)  (7) 

Differentiating  equation  (7)  with  respect  to  time  yields  an  expression 
for  the  pitch  rate, 

A  AQ  /  v 

0  *  Ci>  —  cos  (ojt)  (8) 

From  this  point  the  pitch  rate  can  be  obtained  as  a  function  of  frequency, 
amplitude,  and  neutral  angle  as  follows, 

9  =  2  -  (9  -  Qo)2*  (9) 

The  third  term  of  equation  (6),  g3(h),  represents  the  effect  of  deck, 
height  on  the  induced  pitching  moment  and  has  already  been  discussed. 

For  the  roll  direction,  induced  rolling  moment  can  be  defined  in  a 
similar  manner.  Induced  rolling  moment  is  a  function  of  roll  angle, 
roll  rate,  and  deck  position  as  follows, 


|f  -  H(4>)  +  k2($)  +  k3(h) 


(10) 


Roll  angle  and  rate  are  found  using  the  same  approach  as  was  used  for 
pitch  angle  and  rate.  This  approach  yields  the  following  equations, 


^  sin 


cos  (oot) 


(11) 


(12) 
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The  last  term  of  equation  (12),  k^Ch),  represents  the  effect  of  height 
on  the  induced  rolling  moment  and  has  already  been  discussed. 

DATA  REDUCTION 

This  section  of  the  report  describes  how  the  data  was  reduced  for 
analysis  so  that  a  computer  program  could  be  used  for  computational 
purposes.  The  data  selected  for  analysis,  as  previously  indicated,  were 
obtained  from  reference  (a)  for  the  3-fan  subsonic  VSTOL  configuration. 
These  data  were  in  the  form  of  time  histories  for  induced  lift,  pitching 
moment,  and  rolling  moment.  A  computer  program,  reference  (b),  was  used 
in  order  to  set  up  a  table  of  digitized  values  for  each  data  run.  Once 
established,  this  table  was  entered  and  the  values  at  any  specified  time 
were  determined  using  a  second  computer  program,  reference  (c).  This 
second  program  uses  a  spline  function  to  interpolate  or  extrapolate  the 
table  data  and  has  been  proven  to  be  very  accurate.  These  two  computer 
programs  were  then  used  in  conjunction  with  the  computer  program  REDUCE 
which  was  developed  for  this  analysis. 

DATA  ANALYSIS 

The  computer  program  REDUCE  will  not  be  presented  in  great  detail 
here  but  the  basic  elements  used  to  analyze  the  data  will  be  discussed. 
Figure  A-l  shows  the  flow  of  the  data  analysis.  The  two  points  shown  as 
A  and  B  indicate  the  portion  of  the  flow  for  which  the  computer  program 
REDUCE  was  used.  As  shown  in  this  figure  REDUCE  has  five  basic  elements 
indicated  by  (1)  through  (5).  Element  (1)  is  the  input  of  specified 
rates  for  any  of  the  deck  motions.  These  values  are  determined  for  each 
data  run  based  on  the  equations  previously  developed, 

h/de  *  <u  cos  (tut)  (14) 

T 


Q  ■  u  —  cos  Gut)  (15) 


$  •  a  cos  (ut)  (16) 
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The  maximum  positive  and  negative  rate  values  are  readily  obtained  from 
these  equations  since  the  cosine  has  a  maximum  of  +1  and  the  frequency 
and  amplitudes  are  known. 


Element  (2)  of  this  flow  is  the  computation  of  an  initial  time,  t]_,  which 
is  used  to  compute  all  other  times  of  interest  for  a  given  data  run. 

The  equations  for  rate  were  rearranged  as  follows  to  obtain  t]_. 


t 


1 


t 


1 


1 

—  cos 
0) 


(17) 


(18) 


t.  -  1  cos  4 (19) 

1  u  I  wA<j>| 

The  importance  of  t^  is  that  it  corresponds  to  the  first  time  at  which 
the  specified  rate  occurs  on  the  sine  curve.  As  indicated  by  element  (3), 
all  other  times  at  which  a  specific  rate  occurs  can  be  computed  knowing 
t]_.  The  computation  of  this  value  is  shown  as  element  (4)  of  the  data 
flow  analysis  of  figure  A-l.  Also,  the  number  of  times  that  a  given 
rate  will  occur  in  one  second  must  be  known  and  is  a  function  of  frequency 
as  indicated  in  figure  A-2.  In  this  figure  the  term  hi  is  the  position 
value  (h,  0,  <j>)  corresponding  to  the  specified  rate. 


Element  (5)  is  the  last  segment  of  the  REDUCE  program.  For  each 
time  computed  in  element  (3),  the  induced  force  or  moment  value  was 
obtained  using  the  computer  program  described  in  reference  (c) .  For  the 
entire  5  second  data  run,  an  RMS  of  these  values  was  calculated  along 
with  the  corresponding  position  value.  These  data  were  then  plotted  as 
induced  force  or  moment  versus  position  and  rate  and  have  been  discussed 
in  the  approach  section  of  this  report. 
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Figure  B-5.  Subsonic  Heaving  Motion  Data  (T’vun  364.2) 
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SUBSONIC  PITCHING  MOTION  DATA 
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Figure  C-5.  Subsonic  Pitching  Motion  Data  (Run  82,5) 
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Figure  06.  Subsonic  Pitching  Motion  Data  (Run  82.4) 
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Figure  C-7 .  Subsonic  Pitching  Motion  Data  (Run  85.1) 
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Figure  C-8.  Subsonic  Pitching  Motion  Data  (Run  85.2) 
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Figure  C-9.  Subsonic  Pitching  Motion  Data  (Run  85.3) 
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Figure  C-10.  Subsonic  Pitching  Motion  Data  (Run  87.: 
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Figure  D-l.  Subsonic  Rolling  Motion  Data  (Run  92.1) 
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Figure  D-2.  Subsonic  Rolling  Motion  Data  (Run92.2) 
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Figure  D-5.  Subsonic  Rolling  Motion  Data  (Run  163.2) 


